Self-trapping of Bose-Einstein condensates expanding into shallow optical lattices 
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We observe a sudden breakdown of the transport of a strongly repulsive Bose-Einstein condensate through a 
shallow optical lattice of finite width. We are able to attribute this behavior to the development of a self-trapped 
state by using accurate numerical methods and an analytical description in terms of nonlinear Bloch waves. The 
dependence of the breakdown on the lattice depth and the interaction strength is investigated. We show that it is 
possible to prohibit the self-trapping by applying a constant offset potential to the lattice region. Furthermore, 
we observe the disappearance of the self-trapped state after a finite time as a result of the revived expansion of 
the condensate through the lattice. This revived expansion is due to the finite width of the lattice. 
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I. INTRODUCTION 

The transport properties of Bose-Einstein condensates 
(BECs) through optical lattices have sparked interest in recent 
years after a series of experiments revealed dissipative dynam- 
ics and instabilities [i — f]. Early experiments with BECs in 
optical lattices showed characteristic effects of such a peri- 
odic potential on atoms, namely, Bloch oscillations [^] and 
Josephson junctions [<.,, , ]. In later years the research focus 
shifted towards the study of nonlinear effects arising due to 
interaction of the atoms [S]. Theoretical work suggested that 
two mechanisms, energetic and dynamical instabilities, would 
lead to a dissipative dynamics [9, 10]. This dissipative be- 
havior was also observed experimentally, and it was possible 
to separate dynamical and energetic instabilities [2, 11, 12]. 
In a deep lattice it was found that increasing the nonlinear- 
ity leads to a self-trapped state within the lattice [i . ]. Such 
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FIG. 1: (Color online) Schematic setup of the system. The EEC is 
initially located in the flat-bottom box reservoir A. The shutter to its 
right (dashed vertical line) can be removed instantaneously so that 
the atoms expand into the optical lattice Kpi(z) of size L. A wide 
potential-free region B serves as a second reservoir for the atoms. In 
our numerical calculations we used the dimensionless lengths La = 
IBOtt, Lb = 3267r and L = IOtt (for an explanation of the units see 
text). 
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a self-trapping effect was predicted theoretically in the limit 
of the tight-binding model [14]. A generalized theoretical 
framework for the experiment in Ref. [ 1 3] was derived in [ 1 ?], 
where it was pointed out that the experimental results can be 
explained in terms of a certain type of self-trapped state. 

In this paper, we study the transport of a strongly interact- 
ing, one-dimensional (ID) BEC partially exposed to a shal- 
low optical lattice of finite width. A typical experiment in this 
field is conducted by trapping atoms in a parabolic trap and 
then switching on a moving optical lattice. Alternatively, one 
can displace the parabolic trap in a stationary lattice, which 
leads to acceleration of the atoms through the lattice. We 
model these experiments by a simplified setup where two flat- 
bottom potentials are connected via the lattice as shown in 
Fig. 1. Instead of a moving lattice, we make use of the in- 
herent expansion of the BEC. In contrast to experiments such 
as [l.i], in this present work we focus on the very shallow 
lattice regime, where the lattice depth is strictly less than the 
photon recoil energy. The theoretical study in [ i : ] suggests 
that self-trapping is possible even at lower lattice depths but 
it does not take into account a short extent of the optical lat- 
tice. We extend previous works in this field by assuming an 
optical lattice of finite width connecting two reservoirs. Coun- 
terintuitively, even in the shallow lattice regime, we observe 
self-trapping for sufficiently high interaction strength. In con- 
trast to spatially localized gap solitons [ i : ], this type of lo- 
calization extends over a few lattice sites without decay. The 
localization can be destroyed by applying a constant offset po- 
tential in the lattice region. The numerical analysis also shows 
that the finite width of the lattice results in the dissipation of 
the self-trapping after a finite time. We further compare the 
numerically obtained localized states with an analytical de- 
scription. For our analytical results we incorporate nonlinear 
Bloch waves to approximate the self-trapped state. For the nu- 
merical treatment we simulate the ID Gross-Pitaevskii equa- 
tion (GPE), which is known to be a valid mean-field descrip- 
tion at zero temperature [17]. To study the dynamics of the 
system we use a second-order time-splitting spectral method 
(TSSP) to solve the GPE. This method is explicit, uncondi- 
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tionally stable and spectrally accurate in space. It is known to 
yield higher accuracy for the time evolution of the ID GPE 
with less computational resources compared to, for example, 
Crank-Nicholson methods [ 1 8, 19]. 

One of the underlying motivations of this field is a grow- 
ing interest in the quantum transport of atoms for future mi- 
crodevices. For such a device to become possible it will be 
necessary to identify basic building blocks — the equivalents 
of resistors, capacitors, transistors etc. in electronic devices — 
and to understand their interconnection [20]. One approach is 
to use optical lattices to mimic the crystalline structure found 
in electronic devices [21]. In classical electronic circuits such 
structures are found in wires and many fundamental compo- 
nents such as transistors or diodes. Recently, metallic be- 
havior of ultracold atoms in three-dimensional optical lattices 
has been observed experimentally [ii]. This finding further 
substantiates the analogy between electronic and atomic cir- 
cuits. In an atomic circuit, atomic diodes [23, 24] and tran- 
sistors [ ? ; ] would be connected by atom-transporting wires to 
build more complex devices. Such a wire can be constructed 
with a ID optical lattice, which defines a band structure simi- 
lar to the bands found in metals. In contrast to their electronic 
counterparts, the properties of atomic wires can be dynami- 
cally tuned in a precisely controllable way. By changing the 
angle, frequency or power of the lasers creating the lattice, the 
band structure of the medium can be adjusted to the required 
characteristics. 

In Fig. 1 we sketch our full setup. The flat-bottom reser- 
voirs A and B serve as a source and sink of a BEC, respec- 
tively. The BEC is initially confined to a potential-free box 
A. Such a flat-bottom potential has been implemented experi- 
mentally by Meyrath et al. [26]. The short optical lattice con- 
necting the reservoirs could be implemented by focusing two 
laser beams very tightly. However, we stress that our setup 
is to be understood rather as a generic theoretical model for a 
broader class of experimental setups, for example with shal- 
low harmonic oscillator potentials as reservoirs. The differ- 
ence in chemical potential on the left- and right-hand sides of 
this lattice leads to an expansion of the atoms into the lattice 
and eventually into reservoir B. In the atomic circuit picture, 
the setup can be seen as a battery, and the expansion of the 
atoms leads to a discharge current. 

The paper is organized as follows. In Sec. II we will give a 
detailed overview of the model used and introduce the quan- 
tities of interest. Our numerical and analytical results of the 
transport properties of the BEC are explained in Sec. III. In 
the same section we also briefly discuss the creation of soli- 
tons, which are generated in our scheme. We conclude the 
paper with a summary in Sec. IV. An appendix provides de- 
tails about the calculation of the nonlinear band structure of 
an interacting BEC in an optical lattice. 



II. THE MODEL 

We consider a BEC at zero temperature in the elongated 
trap V{x, y, z) = ^m[Luj_{x^ + y^) + Vax{z)], where m is the 
mass of an atom in the BEC. The transversal frequency u!± is 



chosen to be Hujj_ ^ gno, where g the interaction strength of 
the atoms and uq the peak density of the BEC. In this paper 
we assume 5 > 0, which corresponds to repulsive atomic in- 
teraction. The interaction strength can be expressed in terms 
of the s-wave scattering length as g = AiTfi^as/m. The 
above choice of the frequency uj± results in the freezing of 
the atomic motion in the radial directions. Hence, we treat 
the BEC as an effectively one-dimensional condensate with a 
trapping potential T4x(-z) along the axial direction [27, 28]. In 
our model this potential has the form 




for - La < z < 0, 

for < 2; < L, 

for L < z < L + Lb, 



(1) 



which is iflustrated in Fig. 1. The condensate is initialized 
in the ground state of the box potential of size La, which is 
obtained by solving numerically the time-independent GPE 
using the normalized gradient flow method [29, 30]. After the 
initialization of the BEC in reservoir A, the shutter confining 
the BEC (dashed fine in Fig. 1) is removed. The BEC then 
penetrates a short optical lattice Vopt{z) = Vq + Vi cos{2kz) 
of size L ^ La- We checked numerically that possible dis- 
continuities in T4x at z = and 2 = L do not affect the overall 
results. The periodicity k is given by the geometry and wave 
number of the lasers producing the standing wave, and it de- 
termines the number of lattice sites Lk/ir. The lattice height 
Vi and the constant bias Vb are assumed to be independently 
adjustable. The size of the sink reservoir B is Lb ^ i. 

In order to obtain the dimensionless ID GPE we introduce 
the following dimensionless quantities. Times are rescaled 
according to t = t2E^/h and lengths according to z = kz, 
where L'r = fi^k'^ /2m is the photon recoil energy. Inside 
the lattice region the potential Eq. (1) is given by Wax(z) = 
Vopt{z) = Vopt{z)/ER = V + scos(2z) and zero otherwise. 
Hence, the dimensionless constant offset is v ~ Vo/2Eg^ and 
the lattice depth s = Vi/2Eji. Furthermore, the wave function 
yields 7/>(z, i) = k-'^/'^ipiz, t). At T = the BEC can then 
be described by the dimensionless ID GPE 



idtil' = 



1 



(2) 



where the tildes have been removed for clarity. The dimen- 
sionless interaction strength is /3 = Ntoiaskhuj±/ Er, which 
is expressed in terms of the number of atoms A'tot and the s- 
wave scattering length a^. The wave function tp = -0(2, t) is 
normalized according to J \il-'{z, t)\'^dz = 1 for all times t. 

As an indicator of the dynamics of the system we define the 
dimensionless current 

j{z, t) = 1 [r{z, t)dMz, t) - i,{z, t)d,r {z, t)] . (3) 

As will become apparent in the numerical analysis, it is ad- 
vantageous to also define a more qualitative quantity, namely, 
the stationary current within the lattice. We compute the sta- 
tionary current by taking the time derivative of the particle 
number in reservoir B, A'b, at times where the particle num- 
ber N within the lattice is nearly constant. Given a time 
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FIG. 2: (Color online) Band structure for four sets of parameters ^, s, 
and 71/3. The vertical line separates two sets with the same chemical 
potential (dashed line) but differing s and n/3. Note that each set is 
symmetric around q = 0, hence we plot |g|. The offset v is kept 
at zero. The other parameters are (a) = 0.16 (/3 = 79.58) with 
s = 0.13, n/3 = 0.05 (left part) and s = 0.25, n/3 = 0.04 (right 
part), (b) 11 = 0.795 (f3 = 397.89) with s = 0.095, ?i/3 = 0.329 
(left) and s = 0.127, n/3 = 0.393 (right). The left and right insets 
in (b) show a zoom of the loops near the left and right band edges, 
respectively. 



with such a nearly constant particle number we define the sta- 
tionary cuiTent as 



diVp 



Jo 



dt 



(4) 



to 



In general, the stationary current will depend on all parame- 
ters of the system such as the lattice depth s or the interaction 
strength /3. 



m. BREAKDOWN OF THE ATOMIC EXPANSION 

Ignoring the finite width of the lattice, Eq. (2) reduces to to 
the well-known Mathieu equation in the limit of vanishing in- 
teraction if] = 0). In this limit its energy eigenvalues develop 
the characteristic band structure of periodic potentials. The 
inclusion of the nonlinear term introduces a new energy scale 
n/3 into the system {n is the average particle density within 
the lattice), which changes the band structure. A net effect is 
an overall mean-field shift of the energies by n/3. This effect 
is shown in Fig. 2(a) for typical parameters used in our simu- 
lations. For n/3 > s the band structure additionally develops 
a loop at the band edge, which gradually decreases the width 
of the first band gap [31-33]. These loops can be observed in 
Fig. 2(b). The insets show a closeup of the band edges, where 
a loop has developed. Note that we only plot \q\ in a reduced 
zone scheme, which means that the loop closes symmetrically 
at \q\ > 1. The relative position of the chemical potential 
(dashed line in the figure) and the band gap will be important 
for the explanation of the localized state in the next subsec- 
tion. For details about the calculation of the band structure we 
refer to the appendix. 



A. Numerical results 

In this section we will numerically investigate the transport 
of the BEC initially trapped in region A through the lattice. 



El 











10 150 



10 15 



FIG. 3: (Color online) Upper panel: time-dependent particle num- 
ber within the lattice. Lower panel: time-dependent currents at the 
beginning of the lattice (red solid line) and the last lattice site (blue 
dashed). The gray dashed horizontal line indicates the stationary cur- 
rent obtained as described in the text. Lower panel: The parameters 
are /3 = 397.89, w = for all plots and (a), (c) s = 0.095 and (b), 
(d) s = 0.127. The vertical bar in (b) indicates the analytical result 
for the particle number difference AA*' (see text). 



We will connect the various numerical findings and explain 
them in terms of self-trapped states with finite life time. 

We have calculated the stationary current Eq. (4) numeri- 
cally. Figures 3(a) and (b) show two typical time-dependent 
plots of the particle number within the lattice. We can recog- 
nize plateaux at different times in Figs. 3(a) and 3(b), which 
can be used to compute the stationary current. For example, 
in Fig. 3(b) such plateaux exist at the three time intervals 
[0.008,0.019], [0.021,0.04], and [0.04,0.11]. In the lower 
panel of Fig. 3 two typical results for the time-dependent cur- 
rent Eq. (3) are shown at different positions within the lattice. 
The dashed horizontal line indicates the stationary current cor- 
responding to the same set of parameters. Note that its value 
coincides well with the actual current within the lattice for an 
extended amount of time. The current undergoes small oscil- 
lations around the value of the stationary current. The station- 
ary current indicates the gross expansion speed of the BEC. 
We will analyze its dependence on the parameters s and /3 in 
the following. 

Intuitively one expects the BEC in the setup of Fig. 1 to 
expand into the optical lattice where its transport properties 
are subjected to the modified band structure discussed above. 
We conducted numerical calculations in the regime of weak 
lattices (s <C 1) and strong interaction (/3 ^ 1). If we plot 
the stationary current for different interaction strengths /3 as a 
function of the optical lattice amplitude s, we notice a sharp 
drop in the curves for large f3. In Fig. 4(a) this drop is clearly 
visible, whereas for the lower interactions in Fig. 4(b) it is ab- 
sent in the shallow lattice regime. Note that the value of jo at 
s = increases with increasing /3. This behavior is expected 
because the higher repulsive interaction leads to a higher po- 
tential difference between the reservoirs, which drives more 
atoms through the lattice region. An increase of the lattice 
amplitude does not influence the stationary current in Fig. 4(a) 
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FIG. 4: (Color online) Stationary current for (a), (b) varying s at d = 
and (b), (c) varying constant offset v at fixed s. Parameters for (a) 
are /? = 251.46 (solid line), /3 = 318.31 (dashed), and /? = 397.89 
(dotted), for (b) (5 = 31.83 (solid) and /? = 79.58 (dashed). Plot 
(c) is at /3 = 397.31 for fixed s = 0.095 (sohd line) and s = 0.19 
(dashed), (d) is at /3 = 318.31 with s = 0.095 (sohd) and s = 0.253 
(dashed). The values for the lattice amplitude s used in (c), (d) are 
marked in (a) on the curves with the respective interaction strengths. 



at first, instead it stays constant up to the drop. After the drop 
it decreases with increasing lattice depth. 

We were able to relate the sudden drop in the stationary 
current to the development of an extended plateau in the time- 
dependent particle density within the lattice. In Fig. 3(b) we 
clearly recognize such a plateau for times around 4 x 10^^ to 
11 X 10"^, as well as shorter plateaux at earlier times. For 
the parameters of this plot, (3 ~ 397.89 and s — 0.127, 
a drop in the stationary current has already occurred (cf. 
Fig. 4(a)). Also note that the drop can be observed in the time- 
dependent current in Fig. 3(d), which is lower than the cuiTent 
in Fig. 3(c). This means that after the drop, the BEC den- 
sity within the lattice stays constant for an extended amount of 
time and there is only a small residual current flowing through 
the lattice. The BEC has effectively stopped its expansion 
despite its high repulsion and despite the lattice being very 
shallow. This fact can also be observed in a density plot of 
\'4>{z, t)p. Figure 5(a) shows the density of a BEC with a low 
interaction strength which does not show a drop in the sta- 
tionary current. In contrast, the density for a higher (3 which 
features a drop in the stationary current is plotted in Fig. 5(b). 
The value of s in this plot was chosen larger than the value 
for the drop (cf. Fig. 4(a)). By comparing Fig. 5(b) with the 
particle number in Fig. 3(b) we also identify the steps in the 
plateau structure of the particle number as a tunneling of the 
BEC across lattice sites. Every time a lump of the BEC tun- 
nels to the next lattice site, the average particle density in- 
creases by a fixed amount AA'^ until it reaches a final lattice 
site where the expansion stops and a quasi-stationary state de- 
velops. For a fixed lattice height the final lattice site is de- 
termined by the value of n/3 and s. This becomes clear by 
comparing the chemical potential in Fig. 2 (dashed horizontal 
line) with the position of the first band gap for a given set of 
parameters. For (3 = 79.58, where the stationary current does 
not exhibit a drop, the chemical potential lies deep in the first 
band (cf. Fig. 2(a)). The atoms can populate the first band 
when being injected into the lattice region and move through 
the lattice. For very high (3 however, the chemical potential 
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FIG. 5: (Color online) Density plot of the BEC near the optical lat- 
tice. The lattice extends over < z < IOtt. In the upper panel 
light shades indicate high density, dark shades low density. For a 
typical rubidium BEC such as the one used in [13] the time scale 
of these plots is around 32 ms, which is experimentally accessible. 
The lower panel shows the profile of the BEC at the times indicated 
by the corresponding horizontal lines in the upper panel. The dot- 
ted line indicates the analytical result Eq. (10) for the corresponding 
nonlinear Bloch wave and the thin solid gray line indicates the posi- 
tion of the optical lattice. For both plots s = 0.127, v — and (a) 
f3 = 79.58, (b) 13 = 397.89. 



initially lies above the first band gap (cf. Fig. 2(b)). As we 
increase s, the local chemical potential within the lattice n/3 
also increases slightly and shifts the band structure in such a 
way that the overall chemical potential fj, falls into the first 
band gap. With its chemical potential lying inside the gap 
the wave function cannot expand anymore because there are 
no states available. This manifests itself in the development 
of the quasi-stationary state in Fig. 5(b). Increasing s even 
further broadens the gap and the stationary current decreases 
steadily. We refer to this localized state as self-trapped state 
because its development strongly depends on the value of the 
nonlinearity f3. We find the sudden onset of quasi-stationarity 
within the lattice only for high j3, whereas for lower f3 and 
large lattice amplitude the BEC does not penetrate the lat- 
tice. Alexander et al. [ i ] showed that such a self-trapped 
"gap wave" in a different setup remains stable. However, due 
to the shortness of the lattice in our setup, there is still a small 
cuiTent flowing through the lattice. This leakage eventually 
causes the localized state to disappear after a finite time. In 
Fig. 5(b) this breakdown can be observed for t > 1300, when 
the BEC dissipates over the whole lattice. The phenomenon 
of such a finite life-time of the localized state has also been ob- 
served by Wang et al. in a different optical lattice setup [34]. 

We further studied the case of a fixed s with a varying off- 
set potential v confined to the lattice region. The two curves 
in each of Figs. 4(c) and 4(d) correspond to fixed values of 
(3, respectively. The solid and dashed lines in each plot cor- 
respond to two different values of s. The lattice depths are 
chosen to lie before and after the drop in the stationary cur- 
rent. As expected, for a system whose chemical potential is 
located below the first band gap at v = 0, the stationary cur- 
rent decreases with increasing offset potential (solid lines). In 
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this case the offset shifts the whole band structure by a con- 
stant value. However, if the lattice depth is chosen such that 
the chemical potential lies in a band gap at w = 0, we observe 
a sudden jump in the stationary current (dashed lines). This 
counterintuitive behavior, an increasing current for a higher 
potential barrier, can be explained by again noting that the 
constant offset shifts the band structure. Increasing the offset 
will eventually result in a band structure where the chemical 
potential does not lie in a gap anymore. Thus the jump occurs 
when the chemical potential rejoins a band. 

B. Nonlinear Bloch waves 

The GPE (2) with a periodic potential leads to the well- 
known Mathieu equation in the limit of vanishing interac- 
tion (f3 — 0) if we assume a stationary state -0(2, t) = 
exp(— i/ii)V'q(z). Here, /i is the chemical potential. The 
chemical potential of a BEC in the strongly interacting limit 
can be determined by utilizing the Thomas-Fermi approxima- 
tion, which yields yu = P/La- Our numerical calculations 
of the chemical potential of the initial ID BEC are in good 
agreement with the Thomas-Fermi approximation for the pa- 
rameters used in this paper. The solutions ipq{z) are Bloch 
functions 7/'g(z) = cxp(i(7z) ^^""1"^^ q exp(i2^z). The pa- 
rameter q is the quasi-momentum of the condensate. To model 
the resulting wave function of the interacting case we simi- 
larly assume a Bloch function representation of the state. To 
simplify the analytical model we further truncate the Bloch 
waves according to [33] 

V;,(z) = V^e'-?" (co + + cie^'^) . (5) 

The density n is defined as the averaged relative density of the 
BEC within the lattice 



The normalization of the full wave function 7/)(z) requires that 
|coP + |c_ip+|cip = 1. By using this method it is also possi- 
ble to recover the band structure. For details of the calculation 
of the nonlinear band structure we refer to the appendix. 

To understand the localization of the state after the drop 
consider Eq. (5) with parameters ^/nc_i = \/nci =: di/2 
and do := y/ncQ. For the ground state with q = this results 
in the state 

ip{z) = do + diCos{2z). (7) 
The normalization condition has to be rewritten as 

Mol' + ^MiP-ri. (8) 

The coefficients do and di can be determined by minimizing 
the energy functional under this constraint. This minimization 
procedure yields the two equations 

2f3dl + do{3f3dl - 2n) + sdi = 0, (9a) 
S/3dl + Ml (2 + 3(3dl - /i) + 4s(io = 0, (9b) 



which can be solved analytically. Their real solution describes 
a nonlinear Bloch wave [15, 35]. For the case of g = this 
wave function describes an oscillation with the period of the 
lattice around a finite value. To see this we assume a real 
solution of Eqs. (9) and calculate the density from 
Eq. (7) to yield 

|i/;(z)p = + _L + 2dodi cos(2z) + ^ sin(4z). (10) 

The constant offset in this equation is the particle density n 
given in Eq. (8). The oscillation with the double period 4z can 
be neglected since for our parameter range (i^/2 <C 2|(io'^i|- 
In the simulations we observe that for s above the threshold 
for the drop in the stationary current, a quasi-stationary state 
develops. Figure 5(b) shows the development of such a lo- 
calized nonlinear Bloch wave and its disappearance over time 
for high interaction strength. In the lower panel of Fig. 5(b) 
we show this behavior at two time slices. At early times the 
BEC tunnels across lattice sites as it expands into the lattice. 
At time t w 700 (solid line) we recognize the nonlinear Bloch 
wave in the lower panel of Fig. 5(b). The dotted line overlap- 
ping with the numerical curve is the analytical density Eq. (10) 
with the same parameters as the numerical result. We used the 
reduced chemical potential of the BEC still left in reservoir 
A and the lattice. We note that our analytical description of 
a nonlinear Bloch wave coincides well with the numerical re- 
sult. At a later time t « 1400 (dashed line) the BEC is spread 
uniformly with a periodic modulation. In contrast, for values 
of s and /? where the stationary current does not show a drop, 
the BEC spreads uniformly with time. The lattice region only 
slightly modulates the otherwise uniformly distributed atom 
density. In the lower panel of Fig. 5(a) we see that at t w 30 
(solid line) the BEC has expanded into the optical lattice but 
the density is only slightly modulated with the period of the 
lattice. Similarly, at the later time t « 700 (dashed line) the 
now uniformly spread BEC is only slightly modulated by the 
periodic potential. A quasi-stationary nonlinear Bloch wave 
has not developed. 

We further calculated the size of the steps between the 
plateaux in the particle number plot in Fig. 3(b). Integrating 
the particle density Eq. (8) over one lattice site yields 

A7V = ^(^|do|' + ^Mip) . (11) 

When the BEC advances one lattice site the particle num- 
ber within the lattice should grow by AA^. The vertical bar 
in Fig. 3(b) indicates the analytically calculated difference 
Eq. (11) for the parameters used in the plot with a reduced 
chemical potential as discussed above. This result agrees with 
the difference of the numerically obtained particle number 
plateaux. Given the above agreements of analytical and nu- 
merical results we conclude that indeed the formation of a 
nonlinear Bloch wave causes the breakdown of the stationary 
current. 
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C. Dark solitons 

The GPE supports soliton solutions for nonzero interac- 
tion (3. These solutions are shape-preserving notches or peaks 
in the density which do not disperse over time. In the case 
of repulsive interaction without an optical lattice, solitons are 
typically of the dark type (notches in the density) [36, 37] but 
both dark and bright solitons can exist in BECs in optical lat- 
tices [38, 39]. Our numerical results of the condensate density 
in Fig. 5(b) show the creation of moving dark solitons when 
the condensate jumps to a neighboring lattice site. In Fig. 5(b) 
we can see dark notches moving to the left, away from the 
lattice region. These excitations move slower than the local 
speed of sound (c = \/ph/T^) ^i^^ do not change their shape 
considerably over the simulation time. We also observed other 
typical features of solitons such as the repulsion of two soli- 
tons approaching each other or the phase shift across a soliton. 
Furthermore, we notice that solitons emit sound waves travel- 
ing at the speed of sound. This happens when the center of 
mass of a soliton enters a region of different mean density, 
which causes a change in speed. A detailed analysis of soliton 
trajectories and their deformations in a nonuniform potential 
has been presented in [40]. It should be noted that the creation 
of moving solitons in our simulations shows similarities to the 
creation of solitonlike structures through self-interference of 
the BEC in hard-wall potentials [41]. This self-interference 
could be caused by a small fraction of the BEC being re- 
flected when the larger fraction tunnels through a peak of the 
optical lattice. Furthermore, it is worth noting that the exact 
dynamics and stability of the solitons also depends on the ra- 
tio iJ,/oj± [42]. In our simulations we did not take into account 
the radial confinement, assuming that it is very tight and the 
overall BEC dynamics can be described by a quasi- ID model. 
We did hence not undertake a thorough analysis of the soliton 
dynamics in our system as their creation can be considered a 
side product of the development of the quasi-stationary state 
in the lattice, which was the main focus of this work. 



IV. CONCLUSION 

We have investigated the effects of a finite width lattice on 
the transport properties of a strongly interacting BEC. To this 
end we numerically solved the corresponding ID GPE and 
extracted relevant quantities such as the atomic current and 
density. We also compared the numerical with analytical re- 
sults in terms of nonlinear Bloch waves. 

We found that even for low lattice depths a quasi-stationary 
state develops after an initial expansion of the BEC into the 
lattice. This results in a sharp drop of the current in the lattice 
when the lattice depth and interaction reaches a critical value. 
However, due to the finite extent of the lattice the atoms can 
tunnel out of this state, which eventually leads to the break- 
down of the stationarity. We could explain the development 
of a stationary state with partial nonlinear Bloch waves, which 
builds up over only a few lattice sites and blocks further atom 
flow. When we introduced a constant offset potential into the 
system, we found that increasing the offset can trigger the pre- 



viously suppressed flow of the atoms again. Finally, we re- 
ported on the creation of moving dark solitons during the de- 
velopment of the nonlinear Bloch wave. Every time the BEC 
moves to a neighboring lattice site a soliton is emitted. Hence, 
the number of solitons present in the BEC indicates the num- 
ber of occupied lattice sites. 

In the context of atomic circuits the present work illumi- 
nates the role of wires in an atomic circuit. Our results sug- 
gest that by slightly increasing the lattice depth the superfluid 
current through the wire can be stopped easily. However, this 
effect is not based on a sudden Mott insulator transition since 
we always assumed the applicability of the GPE. In fact, it is 
based on a self-trapped macroscopic configuration, where the 
"charges" penetrate the wire to a certain depth but do not flow 
further. 
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APPENDIX: NONLINEAR BAND STRUCTURE 

As discussed in Sec. IIIB we assume a wave function of the 
form Eq. (5). This trial function can be used to determine the 
energy spectrum of the system. In the noninteracting case this 
leads to the well-known linear band structure. The inclusion 
of the nonlinear term with /? > in the GPE modifies 
the band structure of the Mathieu eigenvalues. The overall 
mean-field shift of the band structure and the development of 
loops have been discussed in Sec. III. Their occurrence can be 
observed by assuming a trial function Eq. (5). We then rewrite 
the normalization condition of the coefficients C£ in terms of 
two angles and 6 according to 

Co — cos 6, (A. la) 

c_i = sin 6' sin 0, (A. lb) 

Cl = sinOcoscj). (A.lc) 

Stationary states can be found by plugging the ansatz 
Eq. (5) together with Eqs. (A.l) into the energy functional 

+ ^i^,(z)r)dz. 

(A.2) 
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This integral is solved analytically and yields lattice potential and the last three lines the interaction energy. 

By fixing different values of the quasi-momentum q and min- 
eg((/), 6) _ 9^ ^ 2 sjj^2 QU ^ cos(2fl!))l imizing Eq. (A. 3) with respect to (f> and 9 we recover the band 

n 2 I- V JJ structure for given parameters n/3 and s (cf. Fig. 2). 

+ - sin(20) [coscj) + sin 0) 

+ — (43-cos(40)[3 + cos(40)] ^^'^^ 
64 I 

-cos(46l)[7 + 8sin(20)] 

+ 8 sin(20) [1 - 008(26*) sin(20)] | . 

The first line represents the kinetic energy, the second line the 
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